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ABSTRACT 

Combustion  instabilities  in  gas  turbine  engines  are  most 
frequently  encountered  during  the  late  phases  of  engine 
development,  at  which  point  they  are  difficult  and  expensive 
to  fix.  The  ability  to  replicate  an  engine-traceable  combustion 
instability  in  a laboratory-scale  experiment  offers  the 
opportunity  to  economically  diagnose  the  problem  more 
completely  (to  determine  the  root  cause),  and  to  investigate 
solutions  to  the  problem,  such  as  active  control.  The 
development  and  validation  of  active  combustion  instability 
control  requires  that  the  causal  dynamic  processes  be 
reproduced  in  experimental  test  facilities  which  can  be  used  as 
a test  bed  for  control  system  evaluation.  This  paper  discusses 
the  process  through  which  a laboratory  -scale  experiment  can 
be  designed  to  replicate  an  instability  observed  in  a 
developmental  engine.  The  scaling  process  used  physically- 
based  analyses  to  preserve  the  relevant  geometric,  acoustic 
and  thermo-fluid  features,  ensuring  that  results  achieved  in  the 
single-nozzle  experiment  will  be  scaleable  to  the  engine. 


INTRODUCTION 

Aggressive  performance  and  emissions  goals  for 
aeroengine  gas  turbine  combustors  have  led  to  the 
development  of  combustor  concepts  that  may  operate  closer  to 


static  and  dynamic  stability  boundaries.  Lean,  direct-injection 
combustors  under  development  share  many  features  with  the 
lean,  premixed  combustors  used  in  stationary  gas  turbines,  in 
which  combustion  instabilities  are  widely  observed. 

In  order  to  minimize  development  costs  and  time,  it  is 
critical  to  posses  the  ability  to  evaluate  the  dynamic  behavior 
of  a combustor  design  during  the  component-development 
phase,  thereby  mitigating  the  need  for  expensive  full-scale 
engine  testing  late  in  the  development  cycle.  Solutions  to 
combustion  instability  problems,  such  as  active  and  passive 
controls,  can  be  evaluated  and  developed  more  quickly  in  a 
laboratory-scale  experiment,  enabling  a faster  transition  of 
technology  into  production  engines. 

Laboratory-  or  component-scale  testing  of  developmental 
combustor  concepts  is  standard  practice  in  the  aircraft  gas- 
turbine  industry,  using  single-nozzle  flame  tubes,  single-  and 
multi-nozzle  sectors  and  full-annular  combustor  rigs.  The  gas 
turbine  community  has  developed  an  experience  base 
regarding  the  fidelity  with  which  such  test  rigs  must  replicate 
engine  design  details  in  order  to  characterize  the  emissions 
and  operability  characteristics  of  engine  combustors.  At  this 
time  there  is  no  proven  methodology  for  replicating  engine- 
scale  combustor  dynamics  in  laboratory-scale  rigs.  The 
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challenge  in  designing  laboratory- scale  combustion  dynamics 
experiments  is  to  replicate  the  engine  dynamic  environment  in 
as  simple  (low-cost)  an  apparatus  as  possible.  Recent 
published  work  suggests  that  bulk  mode  and  longitudinal 
mode  instabilities  can  be  replicated  in  single-nozzle  rigs. 
Cohen  et  al  (1998)  and  Hibshman,  et  al  (1999)  performed 
active  instability  control  experiments  in  single-nozzle  and 
sector  combustors  that  reproduced  a bulk-mode  instability 
observed  in  a lean,  premixed  industrial  combustor.  Paschereit 
et  al  (1998)  have  developed  a sub-scale  combustor  in  which 
the  boundary  conditions  at  the  inlet  and  exit  ends  can  be 
varied  to  impose  a desired  acoustic  mode.  No  relevant  work 
has  been  published  on  replication  of  tangential  modes  in 
multi-nozzle  sector  or  configurations  other  than  full  annula 
rcombustors. 

While  the  aforementioned  investigators  have  shown  that 
it  is  possible  to  create  a realistic,  laboratory-scale  combustion 
dynamics  experiment,  the  scaling  processes  through  which 
these  experiments  can  be  designed  have  not  been  described 
with  sufficient  specificity  to  allow  the  methodology  to  be 
adopted  by  the  technical  community  at  large.  In  order  for  the 
results  of  sub-scale  experiments  to  be  useful  to  engine 
designers,  there  must  be  a system  in  place  to  translate  those 
results  to  engine-scale.  This  paper  discusses  such  a process 
and  presents  an  example  problem,  beginning  with  the  analysis 
of  an  engine-traceable  instability  and  culminating  in  the 
comparison  of  the  lab- scale  results  to  those  from  the  engine. 
The  process  involves  the  following  steps: 

1.  Analysis  of  dynamic  data  from  the  subject  engine  to 

determine  characteristics  of  instability:  frequency, 

amplitude,  and  sensitivity  to  changes  in  hardware 
configuration  and  operating  conditions. 

2.  Acoustic  analysis  of  the  engine  in  order  to  determine 
acoustic  modes  associated  with  instability. 

3.  Conceptual  design  of  single-nozzle  experiment, 
reproducing  the  engine’s  acoustic  environment  and 
replicating  the  relevant  dynamic  processes,  as  determined 
in  steps  1 and  2. 

4.  Acoustic  analysis  of  sub-scale  experiment  to  determine  its 
fundamental  acoustic  modes  in  order  to  confirm  similarity 
with  those  observed  in  the  engine. 

5.  Test  of  the  finalized  laboratory-scale  experiment  and 
comparison  of  data  to  analyses  and  engine  data. 

ENGINE  PROBLEM  ANALYSIS 

Tire  sample  problem  that  was  selected  for  this  effort  was 
a combustion  instability  that  was  observed  during  the 
development  phase  of  a high-performance  aero-engine  that 
employed  a full-annular  combustor  employing  24  fuel  nozzles. 


The  frequency  of  the  instability  varied  from  about  420  Hz  at 
low-power  conditions  to  about  580  Hz  at  high-power 
conditions,  as  shown  in  Fig.  1.  At  a mid-power  operating 
point,  corresponding  to  that  used  for  the  analytical  phase  of 
this  study,  the  frequency  of  the  instability  was  525  Hz.  The 
magnitude  of  the  pressure  oscillations  resulting  from  the 
instability  were  sufficient  to  cause  unacceptable  vibratory 
stresses  in  the  turbine  component  of  the  engine 

Passive  means  of  reducing  the  magnitude  of  the  pressure 
oscillations  were  developed  during  engine  testing.  This  series 
of  tests  was  able  to  detect  a notable  difference  in  the  pressure 
oscillation  magnitude  when  the  fuel  nozzle  /air  s wirier  design 
was  changed. 

A limited  series  of  non-reacting  experiments  were 
conducted  (Anderson,  et  al,  1998)  to  determine  the  acoustic 
and  fluid  mechanic  response  characteristics  of  the  air  swirlers 
that  were  used  in  the  engine  tests.  The  aim  of  these 
experiments  was  to  measure  differences  in  the  responses  of  the 
air  swirler  / fuel  injector  designs  that  might  explain  their 
behavior  in  engine  tests.  It  was  found  that  the  air  swirlers 
associated  with  the  largest  pressure  oscillations  exhibited  a 
preferential  response  to  air  flow  perturbations  in  the  frequency 
range  between  300  Hz  and  500  Hz.  Air  swirlers  that 
performed  better  in  the  engine  did  not  exhibit  this 
phenomenon  in  the  non-reacting  tests. 

While  fast-response  combustor  pressure  data  were 
acquired  during  the  engine  tests,  they  were  acquired  at  a 
limited  number  of  locations.  For  this  reason,  it  was  difficult  to 
draw  any  significant  conclusions  about  the  nature  of  the 
instability  purely  from  the  engine  data.  The  analyses 
described  in  the  next  section  of  this  paper  were  used  to 
augment  and  interpret  the  engine  data.  The  combination  of  the 
engine  data,  the  acoustic  analyses  and  the  non-reacting 
s wirier/injector  characterization  provided  the  basis  for 
replicating  the  problem  in  a single-nozzle  combustor  rig. 


F nsqueiicy,  Hz 
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Figure  1 . "Waterfall"  plot  of  combustor  pressure  spectra 
from  an  engine  test,  showing  evolution  of  the  instability 
during  an  acceleration  event. 

Analysis  of  Engine  Acoustics 

A quasi- ID  unsteady  Euler  analysis  (Wake  et  al.,  1996) 
was  used  to  predict  the  bulk  and  longitudinal  modes  of  the 
engine  combustor.  The  one- dimensional  Euler  equations  were 
solved  with  area  variation  and  modeled  source  terms  to 
account  for  mass  addition,  heat  addition  and  pressure  losses. 
The  resonant  acoustic  frequencies  of  the  combustor  were 
determined  by  solving  for  the  unsteady  time-accurate  response 
of  the  system  and  monitoring  the  fluctuating  pressure  at  a 
specified  location.  The  analysis  can  be  used  to  examine  the 
sensitivity  of  the  frequency  and  pressure  amplitude  to  system 
parameters  such  as  physical  dimensions,  temperature 
distribution,  Mach  number,  flow  rate  and  pressure  drop.  The 
solver  was  first  used  to  compute  the  steady-state  results,  and 
then  unsteady  results  were  obtained  by  forcing  the  system. 
Forcing  was  accomplished  by  adding  an  unsteady  component 
to  the  heat  release.  Two  types  of  forcing  signals  were  used;  1) 
broadband  white-noise  distributed  forcing,  and  2)  swept  sine 
forcing  at  a fixed  frequency.  The  resulting  pressure  response 
indicated  the  frequency  dependence  of  the  combustion  system. 

It  should  be  recognized  that  this  analysis  is  useful  for  its 
intended  role  of  the  determination  of  acoustic  modes.  The 
analysis  does  not  attempt  to  incorporate  the  physics  of  the 
interaction  of  the  fluid  mechanics/acoustics  with  the  heat 
release  (effects  of  flame  shape  variation,  local  fuel/air 
variations,  etc.)  and  therefore  is  not  capable  of  predicting  the 
absolute  magnitude  of  the  pressure  response.  The  approach 
demonstrated  in  this  paper  will  show  that  it  is,  however,  a 
satisfactory  tool  for  evaluating  test  rig  acoustic  characteristics 
for  combustion  dynamics  experiments. 


Figure  2.  Cross-sectional  area  vs.  axial  position  for  quasi-ID 
Euler  model. 

Engine  configuration  Rig  configuration 


Figure  3. Computed  power  spectrum  of  combustor  pressure 
for  engine  configuration  with  sine-sweep  forcing. 

Quasi- ID  Euler  calculations  were  conducted  for  the 
engine  configuration  at  an  intermediate  operating  condition: 
771  F (684K)  and  200  psia  (1.2  MPa).  The  engine  geometry 
was  converted  into  a one-dimensional  description  of  area  vs. 
axial  position  as  shown  by  the  dashed  line  in  Fig.  2,  which 
shows  the  distributions  for  both  the  engine  and  the  rig  (to  be 
discussed  later).  The  geometry  used  included  an  inlet  plenum, 
the  engine  prediffuser,  diffuser  plenum,  the  cowl  or  hood,  the 
swirler,  the  combustor  liner  and  turbine  vanes.  The  combustor 
lies  between  x = 0 and  9.25  in.  Beyond  the  turbine  vane  exit 
the  area  was  expanded  rapidly  to  create  a plenum  dump.  The 
boundary  conditions  used  were  constant  total  pressure  at  the 
inlet  plenum  and  constant  static  pressure  at  the  exit  plenum. 

The  acoustic  response  of  the  system  was  obtained  by 
swept-sine  forcing  of  the  entire  heat  release  distribution.  The 
fluctuation  levels  imposed  on  the  heat  release  were  10%  of  the 
mean.  The  unsteady  pressure  amplitude  at  the  x = 3 in. 
location  in  the  combustor  (3  in.  downstream  of  the  combustor 
dump  plane)  was  used  to  determine  the  pressure  response.  The 
swept-sine  response  is  shown  in  Fig.  3 and  indicates  a 
resonance  at  approximately  575  Hz.  The  width  of  the 
amplitude  response  peak  indicates  a large  amount  of  damping. 
However,  the  width  of  the  response  also  indicates  a broad 
range  of  frequencies  over  which  the  system  may  be 
susceptible  to  combustion  instability. 

These  results  indicated  a longitudinal  mode  in  the 
combustor  near  575  Hz,  near  the  observed  frequency  of  about 
525  Hz.  The  pressure  mode  shape  for  the  550  Hz  mode  is 
shown  in  Fig.  4.  The  mode  represents  a full- wave  solution  to 
the  system  equations  with  zero  unsteady  pressure  specified  at 
each  end.  Given  the  high  impedance  at  each  end  due  to  the 
high  Mach  number  boundaries,  the  mode  shape  can  also  be 
inteipreted  as  a half-wave  across  the  diffuser-combustor 
domain  with  closed  ends.  Note  there  is  a pressure  node 
apparent  at  the  air  swirler/fuel  injector  (x=0).  The  calculated 
fluctuating  pressure  in  the  diffuser  was  1 80  deg.  out  of  phase 
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from  the  pressure  in  the  combustor.  It  is  believed  that  this  was 
the  basic  acoustic  mode  that  occurred  in  the  engine 
configuration  instability.  Although  tangential  acoustic  modes 
exist  in  the  full-annular  engine  combustor,  analysis  of  the 
engine  data  and  two-dimensional  Euler  results  have  indicated 
that  they  are  not  associated  with  the  observed  instability. 


Figure  4.  Computed  pressure  mode  shape  for  575  Hz  mode 
for  the  engine  configuration  at  evaluation  point  conditions. 

SUB-SCALE  EXPERIMENT  DESIGN 

Test  Apparatus  Features 

Based  on  prior  experience  (Peracchio  et  al,  1998)  the  test 
rig  design  approach  incorporated  the  following  guidelines: 

A.  Use  of full-scale  fuel  preparation  sub-components  (fuel 
nozzles,  air  swirlers) 

Prior  evaluations  (Anderson  et  al,  1998)  indicated  that 
the  prototype  injector/swirler  exhibited  an  enhanced 
response  near  the  500  Hz  frequency  observed  in  the 
engine  data  so  it  was  important  to  utilize  that  precise 
design.  While  reduced-size  hardware  may  be  of 
interest  in  order  to  minimize  facility  requirements,  this 
approach  was  not  taken,  due  to  the  introduction  of 
uncertainties  associated  with  reduced-scale  flows. 

B.  Reproduction  of  the  longitudinal  acoustic  behavior 

The  critical  elements  in  this  regard  were  the 
combustor,  the  diffuser,  the  pre-diffuser,  and  the  cowl 
(hood)  which  guided  the  diffuser  air  to  the  fuel  nozzle. 
Table  1 shows  a comparison  of  the  geometric  features 
of  the  annular  engine  burner  with  those  of  the  rig. 
Cross-sectional  areas  were  chosen  to  replicate 
associated  volumes  and  expansion/contraction  ratios. 


Feature 

Engine 

Rig 

Combustor  Volume  per  Injector 
(in:7cc) 

113/ 

1851 

108/ 

1770 

Combustor  Length  (in./eni) 

8.5/21.3 

8.5/21.3 

Shroud  Volume  per  Injector 
(in1  fee) 

129.3/ 

2119 

103.4/ 

1695 

Nominal  Shroud  Height  (in./cm) 

1.2/ 3.0 

0.80/2.0 

Diffuser  length 

2.7/ 6.9 

2.7  / 6.9 

Pre-diffuser  Length  (in./cm) 

3.9  / 9.9 

3.9/ 9.9 

Table  1.  Comparison  Engine  and  rig  acoustic  features. 


C.  Acoustic  isolation  of  the  combustor  from  facility  air 

piping 

A venturi  was  used  to  choke  and  meter  the  inlet  air 
flow.  Since  the  isolation  provided  by  the  sonic  throat 
condition  was  desired  over  a range  of  conditions,  the 
venturi  was  designed  to  be  underexpanded,  resulting  in 
a normal  shock  at  a distance  of  1.38  in.  (3.5  cm) 
downstream  of  the  venturi  throat.  The  upstream 
boundary  was  largely  established  by  the  normal  shock 
and  the  sudden  expansion  of  the  flow  at  the  pre- 
diffuser  dump.  The  downstream  boundary  was  defined 
by  using  a choked  exhaust  nozzle  at  the  station 
occupied  by  the  first  turbine  inlet  vane. 

D.  Reproduction  of  the  airflow  distribution , pressure  drops 

and  flow  damping  characteristics 
The  pressure  drops  and  airflow  splits  used  in  the 
engine  were  duplicated.  That  is,  the  fractions  of  air 
used  for  liner  cooling  and  for  primary  and  dilution  air 
were  reproduced.  Designing  for  equivalent  damping  is 
important  to  achieving  similar  instability  amplitudes 
between  the  test  rig  and  the  engine.  The  resistive 
damping  of  the  system  was  maintained  by  replicating 
the  system’s  pressure  drops.  Reproducing  the  airflow 
splits  between  the  different  air  flows  replicates  the 
distribution  of  stoichiometry  and  heat  release  rate 
within  the  combustor. 

E.  Design  for  testing  at  engine  operating  conditions 

The  instability  observed  in  the  engine  occurred  over  a 
range  of  test  conditions.  A single  “Evaluation  Point” 
was  chosen  at  200-psia  (1.2  MPa)  combustor  pressure, 
an  entrance  temperature  of  771  deg  F (684K),  and 
combustor  fuel-air  ratio  of  approximately  0.03.  All 
analyses  were  conducted  at  these  conditions. 
Operating  at  reduced  conditions  with  full-scale 
hardware  can  change  the  operating  characteristics  of 
the  components  (pressure  drops,  atomization,  etc.). 
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The  final  consideration  was  whether  to  utilize  an 
apparatus  with  a simple  cylindrical  cross-section  burner  or  an 
apparatus  having  a cross-section  representative  of  1/24  of  the 
24-nozzle  engine  burner  - i.e.,  a ‘single  sector’  burner.  A 
circular  cross-section  was  employed,  as  this  represented  the 
lower-cost,  higher- strength  approach.  The  number  and  size  of 
the  combustion  and  dilution  air  holes  was  adjusted  to  provide 
proper  penetration  of  these  air  jets. 


Figure  5.  The  combustor  test  section  assembly. 
Dimensions  are  in  inches. 

The  configuration  of  the  test  section  is  illustrated  in  Fig. 
5.  Provisions  for  high-response  pressure  transducers  and  for 
gas  sampling  (not  reported  herein)  were  incorporated. 

Three  transducer  bosses,  equally  spaced  around  the 
circumference,  were  located  in  the  primary  combustion  zone. 
One  was  located  in  the  secondary  zone,  and  one  in  the  dilution 
zone.  Floating  seals  on  the  liner  captured  extension  tubes 
attached  to  the  pressure-vessel-mounted  transducers  such  that 
the  pressures  at  the  inside  liner  surface  could  be  recorded. 

Bosses  for  shroud  flow  pressure  measurements  were 
provided  at  the  location  of  the  liner  primary  and  dilution  holes. 
Bosses  for  diffuser  pressure  measurements,  upstream  of  the 
combustor,  were  also  provided 

PCB  piezo-electric  pressure  transducers  (P/N  124A21) 
were  selected  for  this  application  based  on  demonstrated 
ability  to  withstand  the  severe  environment  encountered  m the 
combustor  and  pre-diffuser  (high  temperature  and  pressure). 
Originally  developed  for  measuring  combustion  instability  in 
liquid  rockets,  they  are  capable  of  measuring  +250  psi. 
pressure  fluctuations  at  frequencies  between  0.5  Hz  and  10 
kHz.  Satisfactory  durability  was  achieved  by  use  of  an 
integral  water-cooled  mounting  fixture  that  maintained  an 
acceptable  temperature  around  the  sensor.  The  liner  pressure 
sensors  communicated  with  the  combustor  through  a 0.062  in. 
(1.6  mm)-diameter,  0.83  in.(2.1  cm)  -long  sensing  tube.  A 


1000  psia  nitrogen  supply  was  used  to  purge  the  tube  - the 
amount  of  purge  flow  rate  was  negligible.  The  %-wave 
resonant  frequency  of  the  cavity  within  this  tube  was  far  above 
the  frequency  range  of  interest  in  this  experiment.  Analog 
data  were  low-pass  filtered  at  2 kHz  and  digitally  sampled  at  5 
kHz  using  a simultaneous  sample-hold  data  acquisition 
system. 

Analysis  of  Single-Nozzle  Rig  Acoustics 

The  single-nozzle  combustor  rig  design  was  established 
to  preserve  the  axial  lengths  and  cross-sectional  areas  of  the 
engine  configuration  relative  to  a single  nozzle.  The  area  vs. 
axial  position  distribution  was  maintained  approximately  the 
same  (ref.  Fig.  2),  but  some  variation  existed  due  to 
differences  in  engine  hardware  and  the  axisymmetric  hardware 
to  be  used  in  the  single-nozzle  combustor.  Airflow  splits  and 
pressure  losses  - swirler,  bulkhead,  liner,  primary  and  dilution 
jets  were  also  preserved  by  design.  The  inlet  and  exit  of  the  rig 
were  choked  to  acoustically  isolate  the  system.  The  Euler 
models  included  inlet  and  exit  plenums  upstream  and 
downstream  of  the  choke  points  to  allow  constant  total 
pressure  and  constant  static  pressure  to  be  specified, 
respectively,  as  boundary  conditions  to  be  applied  to  the  Euler 
code  domain. 

Quasi- ID  Euler  calculations  were  conducted  for  the 
baseline  rig  configuration  at  the  evaluation  point  operating 
condition:  771  F (684  K),  200  psia  (1.2  MPa).  Swept-sine 
forcing  over  the  frequency  range  from  100  to  800  Hz  was 
applied  to  the  heat  release.  The  resulting  power  spectrum  of 
the  pressure  response  is  shown  in  Fig.  6,  indicating  the 
presence  of  resonances  at  115  Hz  and  550  Hz.  The  level  of 
forcing  employed  in  the  analysis  is  arbitrary,  and  within  the 
linear  response  range,  so  that  no  significance  should  be 
attached  to  the  absolute  levels  of  the  ordinates  in  Figs.  4 and  6. 


Figure  6.  Computed  power  spectrum  of  combustor  pressure 
at  x = 3.0  in.  for  the  rig  configuration.  Quasi-ID  Euler 
code  results  for  100-800  Hz  swept-sine  forcing  of  heat 
release. 
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Figure  7.  Computed  pressure  mode  shape  for  the  550  Hz 
mode.  Rig  configuration  at  evaluation  point  conditions. 

Further  analysis  of  the  mode  shapes  associated  with 
these  resonances  revealed  that  the  low  frequency  1 1 5 Hz  mode 
was  a first-order  longitudinal  mode  where  the  diffuser  and 
combustor  were  in  phase.  The  primary  mode  of  interest  was 
the  550  Hz  mode  since  the  observed  instability  frequency  in 
the  engine  was  525  Hz.  The  pressure  mode  shape  is  shown  in 
Fig.  7.  The  550  Hz  mode  was  essentially  a half-wave 
longitudinal  mode  considering  closed/closed  acoustic 
boundary  conditions  from  diffuser  inlet  to  combustor  exit.  A 
pressure  node  appeared  to  occur  at  the  air  swirler/  fuel  injector 
location.  The  pressure  in  the  diffuser  was  180  deg.  out  of 
phase  from  the  pressure  in  the  combustor.  Note  there  was 
some  activity  downstream  of  the  combustor  exit,  but 
calculations  performed  with  varying  exit  plenum  length  did 
not  indicate  significant  changes  in  the  resonant  frequency. 

The  main  conclusion  drawn  from  the  acoustic  analyses  is 
that  the  single-nozzle  combustor  rig  configuration  will  have  a 
longitudinal  acoustic  resonance  at  about  550  Hz  that  is  very 
similar  to  the  mode  observed  in  engine  data  and  predicted  by 
Euler  analysis  of  the  engine  configuration. 


EXPERIMENTAL  RESULTS 

The  operating  conditions  of  the  combustor  could  be 
completely  described  by  the  following  parameters:  diffuser  air 
pressure  (P3),  diffuser  air  temperature  (T3)  and  combustor 
fuel/air  ratio  (f/a).  Values  for  each  of  these  parameters  were 
chosen  to  correspond  to  three  different  engine  operating 
conditions.  These  are  shown  in  Table  2.  The  fuel/air  ratio 
referred  to  is  that  estimated  at  the  exit  of  the  combustor  and 


accounted  for  all  of  the  air  flowing  into  the  combustor  through 
the  air  swirler,  primary  and  dilution  holes  and  liner/bulkhead 
cooling  passages.  It  was  not  possible  to  vary  the  test 
parameters  independently  because  of  the  choked,  fixed-area 
combustor  exit. 


Inlet  Air 
Pressure,  P3 
(psia/MPa) 

Inlet  Air 
Temperature,  T3 
(F/K) 

Fuel/ Air 
Ratio 

70  / 0.48 

500  / 533 

0.016 

110/0.76 

600  / 589 

0.024 

175/1.21 

771/684 

0.030 

Table  2.  Test  conditions  corresponding  to  engine  operating 
points.  Evaluation-point  conditions  are  in  italics. 

For  the  evaluation  point  operating  conditions,  an 
instability5  was  observed  at  a frequency  of  566  Hz  (see  Fig.  8). 
The  amplitude  of  this  mode  at  these  conditions  was  +/-  0.39 
psi  (2.7  kPa).  The  unsteady  pressure  results  presented  here  are 
from  the  transducer  located  at  1.9  in.  (4.7  cm.)  downstream  of 
the  combustor  bulkhead.  The  amplitude  of  the  instability 
increased  with  increasing  fuel/air  ratio  for  fixed  P3  and  T3. 
At  higher  fiiel/air  ratios,  the  overall  RMS  pressure  fluctuations 
were  dominated  by  this  single  tone.  There  was  significant 
noise  generated  in  the  100-300  Hz  range,  although  none  of  it 
was  particularly  coherent. 


0 500  1000  1500  “ 2000  2500 

Frequency  (Hz) 


Figure  8.  Measured  power  spectrum  of  unsteady  combustor 
pressure  at  x = 1.9  in.  for  Evaluation  Point  operating 

conditions,  showing  resonance  at  566  Hz#  with  an 
amplitude  of  0.39  pi  (0.78  psi  p-p). 
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Figure  9.  Measured  distribution  of  566  Hz  mode,  showing 
magnitude  and  phase  referenced  to  pressure  measurement 
PLA1C1  at  x = 1.9  inches  downstream  of  the  combustor 
bulkhead. 


Figure  10.  Power  spectra  of  unsteady  combustor  pressure 
at  three  operating  conditions  corresponding  to  those  in 
Table  2 showing  decreasing  amplitude  of  566  Hz  mode 
with  decreasing  power  level. 

Figure  9 shows  the  spatial  distribution  of  the  unsteady 
pressure  at  three  locations  within  the  combustor  and  one 
location  in  the  diffuser  region  upstream.  Within  the 
combustor,  there  was  no  phase  difference  between 
measurements  at  different  axial  stations,  and  only  small 
differences  in  amplitude.  Significant  566  Hz  signal  was  also 
apparent  upstream  of  the  combustor,  which  lagged  the 
combustor  pressure  by  92  degrees  in  phase  and  was  smaller  in 
amplitude  by  a factor  of  2.  There  were  no  phase  or  magnitude 
differences  (at  566  Hz)  between  pressure  measurements  at 
equivalent  axial,  but  differing  circumferential  stations, 
consistent  with  a longitudinal  acoustic  mode.  The  unsteady 
pressure  in  the  shroud,  just  outside  the  dilution  holes  was  a 


factor  of  2 smaller  than  that  in  the  combustor  and  lagged  the 
combustor  pressure  by  approximately  40  degrees 

This  mode  was  also  observed  at  the  other  two,  lower- 
power,  operating  conditions,  although  at  a smaller  amplitudes 
and  lower  frequencies,  as  shown  in  Fig.  10. 

ANALYSIS  OF  RESULTS 

Comparison  with  Analytical  Results 

The  experimental  results  can  be  compared  to  the 
analytical  results  by  referencing  Fig.  6,  which  shows  the 
predicted  pressure  spectrum,  and  Fig.  8,  which  shows  the 
measured  pressure  spectrum.  Recall  the  Euler  code  model 
predicted  broad  acoustic  resonances  at  about  575  Hz  and  1 15 
Hz.  In  the  experiment,  a broad  instability  centered  near  570 
Hz  was  observed,  and  there  was  some  incoherent  activity 
indicated  near  1 00-200  Hz.  Thus,  the  agreement  appears  to  be 
good.  The  mode  shape  measured  in  the  experiment  was  of 
limited  spatial  resolution  and  showed  little  spatial  variation  of 
unsteady  pressure  amplitude  or  phase  within  the  combustor 
chamber  itself  for  the  566  Hz  mode  (ref.  Fig.  9).  This  was 
consistent  with  the  mode  shape  of  the  575  Hz  mode  predicted 
by  the  Euler  code  (ref.  Fig.  7).  Note  both  results  did  indicate  a 
slight  decrease  in  amplitude  towards  the  upstream  end  of  the 
combustion  chamber.  The  Euler  code  prediction  indicated  the 
unsteady  pressure  in  the  diffuser  upstream  of  the  combustor 
would  be  180°  out-of-phase  with  the  combustor  pressure.  The 
experimental  results  indicated  a significant  phase  shift  in  the 
diffuser  section,  lagging  the  combustor  pressure  by  about  90 
degrees  at  566  Hz.  This  discrepancy  is  likely  associated  with 
the  1-D  limitations  of  the  model.  For  example,  it  is  expected 
that  some  level  of  coupling  to  the  outer  shroud  passage  would 
occur  which  is  also  coupled  to  the  combustor  via  the  air 
mixing  holes.  Therefore,  some  transition  of  the  phase  from  in- 
phase  with  the  combustor  outside  the  mixing  holes  to  out-of- 
phase in  the  diffuser  section  is  expected  in  the  3-D  problem. 
The  result  could  be  a phase  relation  in  the  diffuser  section 
between  0 and  180  degrees. 

Note  that,  because  the  Euler  code  is  essentially  an 
acoustic  calculation,  it  is  fundamentally  limited  in  its  ability  to 
calculate  the  amplitude  of  the  pressure  oscillations  without  the 
addition  of  a combustion  / acoustic  coupling  model.  In 
calculations  for  the  engine  using  a constant  relative  forcing 
level,  the  Euler  code  indicated  that  both  the  frequency  and 
amplitude  of  the  instability  should  increase  with  increasing 
engine  power  level.  This  trend  was  validated  with  engine  data 
and  was  also  reproduced  in  the  single-nozzle  experiment  (ref. 
Fig.  1 and  Fig.  10).  It  is  also  encouraging  that  the  damping 
mechanisms  present  in  the  calculations  produced  a broad  peak 
at  575  Hz,  much  like  that  seen  in  the  experiment  (see  Fig.  10). 
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Comparison  with  Engine  Data 

Figure  11  shows  a comparison  between  the  fluctuating 
pressure  spectrum  in  the  engine  and  the  single-nozzle 
combustor  at  comparable  operating  conditions.  Both  data  sets 
were  acquired  over  10  seconds,  and  were  processed  using  the 
same  techniques.  The  frequency  of  the  “target”  mode  was 
reproduced  within  12%.  The  amplitude  of  this  mode  was 
matched  within  3%.  The  spectral  peak  was  significantly 
narrower  in  the  engine  data,  indicating  a more  coherent 
instability.  The  single-nozzle  combustor  also  exhibited  a 
higher  overall  level  of  noise  in  the  signal,  especially  at 
frequencies  below  350  Hz. 


NASA  SNR  R004 p38  and  Engine  TX3081.2B 


Frequency  (Hz) 

Figure  1 1 . Comparison  of  engine  and  combustor  rig 
pressure  spectra  for  Evaluation-Point  operation. 

CONCLUSIONS  AND  RECOMMENDATIONS 

A methodology  for  replicating  longitudinal  combustion 
instabilities  observed  in  an  aircraft  engine  in  a single-nozzle 
test  rig  was  successfully  demonstrated.  The  experiment 
reproduced  the  frequency  of  the  engine  instability  within  12% 
and  its  amplitude  within  3%. 

Necessarily,  the  replication  procedure  must  use  both 
predictions  from  analytical  tools  and  engineering  judgements 
based  on  prior  combustion  dynamics  studies.  A relatively 
simple,  quasi-one -dimensional  Euler  analysis  is  satisfactory 
for  the  prediction  of  the  dynamics  of  longitudinal  instabilities 
in  geometrically  complex  burners. 

The  one- dimensional  analytical  acoustic  tools  applied  in 
this  program  predicted  the  basic  acoustic  frequencies  of  both 
the  engine  and  the  rig  writhin  about  10%,  and  supported  a 
methodology  to  define  a test  rig  with  specific  acoustic 
frequency  characteristics.  However,  because  the  tool  is  a one- 
dimensional analysis,  it  cannot  capture  all  the  dynamic 


features  of  this  combustor.  For  example,  as  reported,  it  did 
not  quantitatively  predict  the  phase  relationship  between  the 
combustor  and  diffuser.  The  development  of  an  analytical 
tool  wrould  include  multi-dimensional  capability  to  either 
identify  other  modes  or  include  the  influence  of  parallel 
acoustic  paths  (e.g.,  shroud)  is  recommended. 

An  extension  of  these  tools  to  design  sub-scale 
experiments  in  which  instabilities  associated  with  tangential 
acoustic  modes  would  also  be  useful.  Inclusion  of  acoustic- 
heat  release  coupling  would  also  increase  the  tooFs  utility. 
Acoustic-heat  release  coupling  models  could  provide  a 
capability  for  further  reductions  in  the  scale  of  the  experiment, 
in  such  parameters  as  geometric  size,  flow  rates  and  operating 
pressure. 

In  summary,  a methodology  for  design  of  laboratory 
combustors,  which  reproduce  combustion  instabilities 
observed  in  aircraft  engine  combustors,  has  been 

demonstrated.  A realistic  platform  for  the  development  and 
validation  of  active  combustion  instability  control  systems  was 
constructed,  and  will  be  utilized  in  subsequent  programs. 
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PAPER  -5,  W.  Proscia 

Question  (A.  Annaswamy,  USA) 

Is  peak-splitting  due  to  active  control  (in  general)  or  phase-shift  control  (in  particular)? 
Reply 

The  plant  delay  in  the  sector  combustor  is  much  longer  than  would  be  necessary  on  an 
engine,  due  to  long  fuel  valve  to  injector  length.  This  large  delay  made  broad  band  at- 
tenuation of  pressure  oscillations  difficult.  Because  of  the  long  delay  between  the  valve 
command  signal  and  pressure  pulsations,  the  plant  phase  changed  rapidly  over  a rela- 
tively wide  range  of  frequencies  near  the  resonant  frequency,  where  the  plant  had  consid- 
erable gain.  The  quickly  changing  phase  resulted  in  positive  feedback  control  on  both 
sides  of  the  attenuation  band  centered  at  the  resonant  frequency,  yielding  an  apparent 
splitting  of  the  dominant  instability  mode. 

More  information  is  available  in  ASME  Paper  99-GT-215  by  Hibshman  et  al.,  presented 
at  the  ASME  Turbo  Expo  in  Indianapolis,  Indiana,  June  1999. 

Question  (F.  E.  C.  Culick,  USA) 

What  damping  do  you  have  in  the  Euler  code  and  why  do  you  seem  to  have  more  noise 
predicted  by  the  Euler  code  than  measured  in  tests? 

Reply 

The  Euler  code  includes  lumped  losses  to  represent  pressure  losses  in  the  air  swirler  and 
at  the  dump  plane.  Air  flow  through  linear  cooling  holes,  primary  holes,  and  secondary 
holes  is  modeled  as  a plenum  fed  with  appropriate  pressure  losses  across  the  linear. 

Regarding  the  differences  between  the  prediction  and  the  tests,  the  swept-sine  spectra  are 
obtained  by  forcing  the  heat  release  in  the  Euler  code  (Figures  3 and  6 in  the  paper)  and 
not  by  an  unstable  feedback  system.  The  experimental  amplitudes  are  likely  due  to  the 
nonlineaityr  of  the  heat  release  response  leading  to  limit  cycle  behavior.  Also,  the  Euler 
code  spectra  are  plotted  on  a log  scale,  the  experimental  results  on  a linear  scale. 
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